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Abstract 

We adapt White's density matrix renormalisation group (DMRG) to the direct study of critical 
phenomena. We use the DMRG to generate transformations in the space of couphng constants. 
We postulate that a study of density matrix eigenvalues leads to a natural identification of rcnor- 
malised blocks, operators and Hamiltonians. We apply the scheme to the phase transition in the 
anisotropic spin-1/2 Heisenberg chain. In the simplest case where the two most probable states 
in odd sized blocks are used to construct approximate renormalisation group transformations, we 
find qualitative improvement upon the standard real space renormalisation group method for the 
thermal exponent v. 



1 Introduction 

The advent of White's density matrix renormaUsation group (DMRG) method [Q has lead to some 
very successful studies of low-lying excitations and static correlation functions in a number of one 
dimensional quantum lattice systems The method represents a major improvement upon its 
precursor the (conventional, Wilson) real space renormalisation group (RSRG) method [| which 
generally gives poor or slowly converging results for these quantities. 

Both methods are truncated basis expansions in that a target state (such as the ground state) of a 
large lattice is built up from blocks of sites from which only a few important states are retained. The 
key difference between the two methods is the way in which the important states are determined. For 
the RSRG states are retained on the basis of energy whereas in the DMRG states are retained on the 
basis of how likely they are to be part of the target state being investigated. 

It is well documented that the DMRG method works best when the system being studied is away 
from criticality — that is, when there is a substantial energy gap, or when the generic correlation 
functions decay exponentially, with a correlation length of only a few lattice spacings 0, @, ||. In fact, 
efforts to study zero temperature phase transitions by using obvious approaches such as investigating 
the divergence of the correlation length jj] or the vanishing of an order parameter have failed 
because it is very difficult to determine these quantities accurately near the critical point. 

It is well known, however, that the RSRG can be used to generate renormalisation group trans- 
formations in the space of coupling constants Useful qualitative and even accurate quantitative 
results can be obtained for the positions of phase transitions and critical exponents Q . 

A next step is to investigate as to whether the DMRG can be used to generate such transforma- 
tions and yield accurate phase diagrams and associated critical exponents, especially in cases where 
the RSRG fails. That is, it would be useful to develop a DMRG scheme for coupling constant transfor- 
mations which has the same accuracy, systematic improvability and portability as the DMRG scheme 
currently being used for the calculation of energies and correlation functions. Such a scheme might 
be very useful if applied to the investigation of critical phenomena in complex models such as coupled 
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fermion chains. 

Such an investigation has been carried out in a recent series of papers by Drzewihski and coworkers 
Results from the RSRG were compared directly with those obtained from its DMRG analogue. In 
studies of anisotropic XY models with transverse fields, it was generally concluded that the DMRG 
has no special advantage over the RSRG in calculating critical points and exponents. These studies 
cannot however be considered exhaustive. Firstly, the blocks and superblocks used were small — only 
a few sites. Secondly, the models did not afford the total z-spin as a good quantum number. This 
sometimes complicated the process of identifying important block states with renormalised block spin 
variables. 

In this paper we consider a DMRG approach to coupling constant transformations which makes use 
of the fact that within DMRG algorithms we accurately calculate important states for very large blocks 
by retaining large numbers of states at each iteration. We apply this approach to a simple model, the 
anisotropic Heisenberg model, where critical properties are very well known [ p^ . The model has been 
studied using the RSRG with less than encouraging results. We find that the DMRG approach 
yields qualitative improvement in the nature of the convergence of the thermal critical exponent ly as 
the blocksize is increased. 

In the following sections we briefly outline the DMRG algorithm and describe how it can be 
used to generate coupling constant transformations. We then present our results for the anisotropic 
Heisenberg model and compare them with the RSRG results of reference |^ . We then conclude with 
some remarks on future directions. 

2 The DMRG and coupling constant transformations 

The DMRG algorithm has been described in great detail jl] so we will be brief in our description, 
concentrating on those points which are relevant to generating renormalisation group transformations. 
We restrict ourselves to the infinite lattice algorithm which we use in our calculations and we describe 
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the algorithm in the context of the Heisenberg spin chain 

where Si is the spin operator of spin S for site i on the chain and 7 is the anisotropy. 

The DMRG is an iterative, truncated basis procedure whereby a large chain (or superblock) is 
built up from a small number of sites by adding a small number of sites at a time. At each stage 
the superblock consists of system and environment blocks (determined from previous iterations) in 
addition to a small number of extra sites. Also determined from previous iterations are the matrix 
elements of the block Hamiltonians and the active spin operators (those on the sites at the end(s) of 
the blocks) with respect to a truncated basis. Tensor products of the states of the system block, the 
environment block and the extra sites are then formed to provide a truncated basis for the superblock. 
The ground state \ip) of the superblock is determined by a sparse matrix diagonalisation algorithm. 

Next, a basis for an augmented block, consisting of the system block and a specified choice of the 
extra sites, is formed from tensor products of system block and site states. The augmented block 
becomes the system block in the next iteration. However, in order to keep the size of the superblock 
basis from growing, the basis for the augmented block is truncated. We form a density matrix by 
projecting onto the augmented block which we diagonalise with a dense matrix routine. We 

retain the most probable eigenstates (those with the largest eigenvalues) of the density matrix in order 
to form a truncated basis for the augmented block that is around the same size as the system block 
basis. Matrix elements for the Hamiltonian and active site operators are then updated. 

The environment block used for the next iteration is usually chosen to be a reflected version of 
the system block. The initial system and environment blocks are chosen to be a small number of 
sites, usually one or two. The accuracy and computer requirements of the scheme is fixed by ng, the 
number of states retained per block (of good quantum numbers) at each iteration, rig determines the 
truncation error, which is the sum of the eigenvalues of the density matrix corresponding to states 
which are shed in the truncation process. The error in quantities such as the ground state energy 
scale linearly with the truncation error 0] . 
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Now, at each iteration we generate coupling constant transformations T as follows. After forming 
the system block of size L we construct a lattice C of size 2L consisting of two system blocks. A 
(small) subset B of the system block truncated basis states is identified with a complete basis for a 
renormalised block of spins of size L' where L' = L/b and b is the renormalisation factor. B(E)B is then 
clearly identifiable with a complete basis for a lattice C of size 2L'. Next, the matrix elements of H, 
the Hamiltonian for C, with respect to B®B (which are readily formed from the matrix elements of the 
block Hamiltonian and active site operator(s)) are identified with those of a renormalised Hamiltonian 
Ti! defined on C . T is then defined by 

i = r(6|7) (2) 

where 7 and 7' are set(s) of coupling constant (s) which define the Hamiltonians Ti and H! respectively. 

To complete the prescription we must have a suitable method for choosing the set B — it's size and 
makeup — and hence the renormalisation factor 5, together with the type of renormalised block it is 
identified as a basis for. We must then check that our identification is valid in that the Hamiltonian 
matrix element identification can in fact be made such that W and 7i, in a loose sense, have the same 
symmetries. Ideally, we would like 7i' and TL to be of the same form. The resulting system of linear 
equations for the renormalised couplings 7' is typically overdetermined. This implies a consistency 
check of the identification. 

As to the choice of -B, just as the density matrix eigenvalues determine which states are to be 
retained in forming a new system block, so too can we use the spectrum of the density matrix in order 
to guide our choice as to the makeup of B. That is, we choose B from states that make up the bulk 
of the weight of the density matrix whose eigenvalues sum to unity. 

Choice of the basis B 

We commence with the simple example of an S' = 1/2 model with 7=1 (isotropic case). We use 
a periodic superblock of the form . . . -site-system-site-environment-. . . and augment the system block 
with both sites. The strongest density matrix eigenvalues are listed in table 0for the case of ng = 20. 

4 



The initial system block is a single site so the blocks always have an odd number of sites. The 
strongest density matrix eigenvalues lie in the sector of the Hilbert space where the good quantum 
number = Sf has small magnitude. We see that there is a pair of strong eigenvalues in the 
— sectors which make up the bulk of the weight of the density matrix. 

An obvious choice for B then is the states corresponding to these eigenvalues. We identify them 
with the up and down z-spin states of a renormalised spin-1/2 operator. That is, the system block is 
renormalised to a single site — b = L. As we shall see, the advantage of such a choice of B is that Ti' 
is of the same form as Ti. 

Next, in table |^ we list the strongest density matrix eigenvalues for the case of an open ended 
superblock of the form system-site-site-environment. We again augment the system with both sites. 
The initial system block is a pair of sites so the blocks are always even. We see that the simplest 
possible identification involves the four strongest eigenvalues — a pair in the 5^ = ±1 sectors and a pair 
in the = sector. It is natural to form B from the four states corresponding to these eigenvalues 
and to identify _B as a basis for a block of two spin-1/2 sites ie b — L/2. Following we identify the 
elements of B with singlet and triplet states. That is, we identify the — ±1 states with | |t) and 
11) and the strong and weak eigenvalue (low and high energy) — states with -^(| |J.) — | jf)) 
and -^(1 II) -I- 1 ID) respectively. It is found that with this identification 7i' has the same form as H. 

Finally, in table ^ we list the strongest density matrix eigenvalues for a spin-1 model using odd 
sized blocks, periodic boundary conditions, augmenting two sites at a time. We see that, for small L, 
the bulk of the density matrix made up from the three strongest eigenvalues which lie in the = 
and — ±1 sectors. When the corresponding states are (naturally) identified with the |0) and | ± 1) 
states of a spin-1 operator, the resulting H' has the same form as H. However, as L is increased, 
the second largest eigenvalue in the S*^ = sector rapidly becomes comparable to the other retained 
states' and the identification of the renormalised block with a spin-1 site becomes inconsistent. It 
becomes necessary to incorporate the corresponding state into B which, as in the previous example, 
is identified as a basis for a block of two spin-1/2 sites. 
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The renormalised 5=1/2 Hamiltonian TL' is not however of the form ([l|) in that there are second 
and third neighbour interactions present with dimerisation. This is consistent with White's analytic 
mappings between the spin-1 chain and coupled spin-1/2 chains [ p2| . That is, it becomes natural to 
identify the renormalised block with a single segment of couple chains. 

2.1 Discussion — systematic improvability of transformations 

Now that we have described the procedure for generating coupling constant transformations, we discuss 
the question of its accuracy and systematic improvability. These issues are central to the success of 
the DMRG in calculating energies and correlation functions. We consider how the algorithm should 
be scaled in order to obtain the most accurate results. 

Now, suppose we choose B in the same way at each iteration eg for a spin-1/2 system with odd 
blocks we always form B from two states — the most probable states in the — sectors. We have 
the guiding principle that, as long as we respect the basic symmetries of the Hamiltonian, then for 
a fixed basis size \B\, the accuracy of the transformations should increase with each iteration, with 
calculated quantities converging to their exact values as the blocksize is increased. 

This principle stems from the fact that the ratio of intrablock to interblock components of the 
Hamiltonian decreases as the blocksize is increased and the exact ground state can asymptotically be 
written as a suitably symmetrised product of the B states. The principle in fact appears to be borne 
out in the case of the spin-1/2 transverse Ising (ITF) model where slow but systematic convergence 
of critical exponents with blocksize is achieved within a RSRG scheme whereby blocks of sizes 3,5, 
7, 9 and 11 are diagonalised and the basis B is formed from the states of lowest energy 

This principle as applied to estimates of critical exponents is however non-rigorous in general and 
within the DMRG scheme there is likely to be a limit to accuracy of infinite lattice results imposed 
by the finiteness of rig. That is, the finiteness of rig imposes restrictions on how accurately we can 
determine important states of large blocks. In practice it may therefore be profitable to enlarge B 
and hence L' . This may have the undesirable side effect of TC lying in an enlarged space with more 
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interactions. 

Finally, we describe how the approach considered here differs from that used in In rela- 
tively small superblocks (up to 12 sites) are diagonalised exactly, reduced density matrices are then 
constructed for small system blocks (either 3 or 4 sites) and the basis B is formed from retaining the 
2(4) most probable states for 3(4) site system blocks. This is the direct DMRG analogue of the RSRG 
calculations Q. In the approach considered here we again use only a handful of states to generate 
the RG transformation but many states are retained at each iteration for the purpose of constructing 
successive system blocks. That way we attempt to work with accurately determined important states 
of large blocks. Also, the density matrix eigenvalues are used in order to obtain a natural selection of 
the basis B. 



3 Application to the spin- 1/2 chain with anisotropy 

We now apply our method to the 5* = 1/2 case of the Hamiltonian given in equation ([T]) . This model 
is integrable and many properties have been calculated exactly . We only consider 7 > here. In 
this regime the model has two trivial, stable fixed points. At 7 = (Ising fixed point) the ground 
states are the classical Neel states and there is a finite energy gap and long range antiferromagnetic 
order. At 7 = 00 {XY fixed point) the model is equivalent to a spinless fermion gas the spectrum 
is a gapless continuum and there is no long range order. There is a phase transition at 7 = 7c = 1 
(isotropic, Heisenberg point). This critical point separates the gapless phase from the gapped, ordered 
doublet phase. 

The phase transition is pathological in that the critical exponents are either zero or infinite. For 
example, the thermal exponent describing the divergence of the correlation length ^ is usually 
defined by 

17-7^1^" as 7^ 7c (3) 
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However, for the 5=1/2 model we have the exact result Q 

^ ^ ^~A\,-M-''- (4) 

whence v = oo. 

As mentioned, the RSRG method has been applied to the S = 1/2 model |^. Blocks of size L — i, 
5, and 7 were used and the basis B was comprised of the two low energy states (in the S'f = ±1/2 
sectors) and was identified with the up and down states of a spin-1/2 site. The estimates of v so 
obtained were spurious in that v [2 as L — > oo. 

Here we calculate v from RG transformations derived within the DMRG scheme described above. 
We use the periodic superblock, augmenting two sites per iteration with a single site as an initial 
system block (hence L is always odd). We choose = 2 (L' = 1, 6 = L), identifying the most 
probable states in the S'f = ±1/2 sectors, which we denote by | ])) and | |)), with the up and down 
states of a single spin-1/2 site. The renormalised Hamiltonian then has matrix elements 

{a'p'\n'\ap) ^ {{p'\ ® {{cj'\n\a)) ® \p)) (5) 

TL' has sixteen matrix elements \e B ® B has four elements. However, the DMRG algorithm 
preserves certain symmetries of the Hamiltonian at every iteration. From conservation of total z-spin 
we have 

(ctVIHVp) = if o-' + pV + P (6) 
From spin-flip symmetries we also have the following obvious relations 

(TTIH'ITT) - (iilH'lii) (7) 
(TilH'ITi) = (iTiH'UT) (8) 
(TilH'UT) = (iTlH'ITi) (9) 

These symmetries reduce the number of independent matrix elements to 3 ie (tt |^'| TT)i (it I'^'l it) 
and (Ti \n'\ it). 
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The identification of TL' is simple. We can write TL' in the same form as TL (up to a scale factor J' 
an additive constant C") viz 

W = J' [SlS'i + 7' (S-i^S-f + S\Sl)] + C7 (10) 

where / denotes the identity operator and J', C and 7' are uniquely determined from the three 
independent matrix elements. As mentioned, in general the identification is overdetermined. 

It is found that the renormalisation group transformation 7' = T(6|7) converges quite rapidly and 
uniformly with ns . We plot the transformation for various block sizes b (or L) in Fig. |l] in the case of 
ns = 50. We also include a plot of the RSRG result |9| for & = 3 viz 

T(3|7) = (11) 

1 + VI + 87^ 

We see that the RSRG and DMRG results are almost identical for 6 = 3. We shall see, however, that 
a marked difference that emerges for larger h. 

We see that the transformation has trivial fixed points at 7 = and 7 = 00 and the critical point 
7 = 7c = 1 is recovered as an unstable fixed point. This is to be expected, again from symmetry 
considerations ie TL only has full rotational symmetry at 7 = 7c. 

Note that the same argument applies in the spin-1 case. As mentioned, however, it was found that 
the self-mapping for the 5 = 1 case rapidly became inconsistent as b was increased. This is consistent 
with exact diagonalisation results jl^ where it is found that the S — 1 model is not critical at 7 = 1 
but rather at 7 = 1.167(7). 

Results for the thermal exponent p 

Following we use the following standard relation in order to obtain the thermal exponent 

" logr'(6|7c) ^ ' 

We plot the estimates of the exponent as a function blocksize b in Fig. ^ for various values of Tig- 
We also include the RSRG results from We see that the DMRG results are qualitatively more 
consistent with the exact result = 00 in that DMRG result increases with b whereas the RSRG result 



monotonically decreases with b. Further, the DMRG resuh may increase without bound as 6 — ^ oo. 
We plot versus b for the Ug = 50 case in Fig. This function appears to be hnear for large 5, 
consistent with v diverging logarithmically with b. 

4 Discussion 

In this paper we have presented a scheme for using the DMRG algorithm to generate KG transforma- 
tions in the space of coupling constants. The scheme allows us to study critical phenomena directly 
using the DMRG method. The scheme makes use of the fact that the important states of large blocks 
are accurately calculated within the DMRG algorithm. Also, a natural solution to the perennial 
problem of the choice of renormalised block and the identification of the renormalised Hamiltonian is 
proposed in terms of the density matrix spectrum. This allows us to distinguish the different univer- 
sality classes of the isotropic point in the anisotropic spin-S* Heisenberg model. For 5* = 1/2, the RG 
transformation is a self map with the isotropic point as a critical point. For S — I the model maps 
onto an S* = 1/2 model with dimerisation and longer ranged interactions that can be identified with 
a coupled S — 1/2 chain. 

We have applied the method to the pathological phase transition in the 5=1/2 model. Results 
for the thermal exponent are qualitatively better than the comparable RSRG results. That is, the 
exact result is — (x. The RSRG yields | 2 as 5 — > oo (where b is the blocksize) whereas the DMRG 
has V increasing with b, seemingly without bound. 

We believe that the scheme described here may be useful in studying critical phenomena in quasi- 
one dimensional systems with many states per unit cell such as coupled spin or fermion chains, ex- 
tended Hubbard-type models with a charge transfer gap and spin or fermion models with dimerisation 
and frustration. 
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Table Captions 

1. Strongest density matrix eigenvalues for the isotropic 5=1/2 model using rig = 20 for various odd 
blocksizes b using periodic boundary conditions and augmenting by 2 sites per iteration. 

2. Strongest density matrix eigenvalues for the isotropic S = 1/2 model using rig = 20 for various 
even blocksizes b using open boundary conditions and augmenting by 2 sites per iteration. 

3. Strongest density matrix eigenvalues for the isotropic S = 1 model using Ug = 20 for various odd 
blocksizes b using periodic boundary conditions and augmenting by 2 sites per iteration. 



13 



b 




Strongest Eigenvalues 


3 


±3/2 


0000 


3 


±1/2 


0.5000, 0.0000, 0.0000, . . . 


c; 

o 


— 1— o / ^ 




5 


±1/2 


0.4619, 0.0271, 0.0055, . . . 


9 


±:-!/2 


0.0120. . . . 


9 


±1/2 


0.4251, 0.0484, 0.0120, . . . 


21 


±3/2 


0.0237, . . . 


21 


±1/2 


0.3786, 0.0682, 0.0237, . . . 



Table 1: 



14 



b 




Strongest Eigenvalues 


4 


±1 


0.02752, 0.00000, . . . 


4 





0.91746, 0.02752, 0.00000 . . . 


6 


±1 


0.03917, 0.00003, ... 


6 





0.88237. 0.03917. 0.00003 


10 


±1 


0.05204, 0.00022, ... 


10 





0.84310, 0.05204, 0.00022, . . . 


22 


±1 


0.06938, 0.00104, . . . 


22 





0.78808, 0.06938, 0.00104, . . . 


Table 2: 


b 




Strongest Eigenvalues 


3 


±2 


0.00000, . . . 


3 


±1 


0.33333, 0.00000, . . . 


3 





0.33333, 0.00000, 0.00000, . . . 


5 


±2 


0.00520, . . . 


5 


±1 


0.29495, 0.00520, . . . 


5 





0.29495, 0.07181, 0.00520, . . . 


9 


±2 


0.00579, . . . 


9 


±1 


0.26230, 0.00579, . . . 


9 





0.26230, 0.15410, 0.00579, ... 


21 


±2 


0.0038(8. . . . 


21 


±1 


0.23786, 0.00389, . . . 


21 





0.23786, 0.22696, 0.00389, . . . 



Table 3: 
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Figure Captions 

1. The renormalisation group transformation 7' — T{b\'j) for the anisotropic spin-1/2 chain from the 
DMRG method with Ua = 50 for various blocksizes b. The thick dashed curve is the RSRG result for 
b = 3. 

2. DMRG estimates of the critical exponent u as a, function of blocksize b for ris = 20 (dashed Une), 
rig = 35 (dot-dashed hne) and rig = 50 (sohd line). The thick dashed curve is the RSRG result. 

3. e'' versus b for the rig = 50 case. 
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Figure 1: 



Figure 2: 



Figure 3: 
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